function a=difdiv(x,y)

	%functie care calculeaza diferentele divizate pentru o functie data prin valori pe un interval
	% Intrari:
	%  x = abscisele celor n puncte
	%  y = ordonatele celor n puncte
	% Iesiri:
	%  y =diferente divizate de ordin 1:n

	n = length( x );
	for k = 1 : n-1
		y( k+1 : n ) = ( y( k+1 : n ) - y( k ) ) ./ ( x( k+1 : n ) - x( k ) );
	endfor
	a = y( : );
endfunction
